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The mechanism describing the recently developed notion of kernel gravity waves (KGWs) is re- 
viewed and such structures are employed to interpret the unstable dynamics of an example stratified 
plane parallel shear flow. This flow has constant vertical shear, is infinite in the vertical extent, and 
characterized by two density jumps of equal magnitude each decreasing successively with height, 
in which the jumps are located symmetrically away from the midplane of the system. We find 
that for a suitably defined bulk-Richardson number there exists a band of horizontal wavenumbers 
which exhibits normal-mode instability. The instability mechanism closely parallels the mechanism 
responsible for the instability seen in the problem of counter-propagating Rossby waves. In this 
problem the instability arises out of the interaction of counter-propagating gravity waves. We argue 
that the instability meets the Hayashi- Young criterion for wave instability. We also argue that the 
instability is the simplest one that can arise in a stratified atmosphere with constant shear flow. 
The counter propagating gravity waves mechanism detailed here explains why the Rayleigh criteria 
for shear flow instability in the unstratified case does not need to be satisfied in the stratified case. 
This illustrates how the Miles-Howard theorem may support destabilization through stratification. 
A normal mode analysis of a foamy layer consisting of two density jumps of unequal magnitude is 
also analyzed. The results are considered in terms of observations made of sea-hurricane interfaces. 



I. INTRODUCTION 

Shear flows and the variety of instabilities they precipi- 
tate are difficult processes to conceptualize despite over a 
century of inquiry. The two different interpretative tools 
available are the theory of counter-propagating Rossby 
Waves (CRWs)i2j2i and classic over-reflection theory (O- 
R)^ and these, in turn, have been shown to be equiv- 
alent in rationalizing non-stratified shear flows 5 . While 
O-R theory can be used to interpret stratified shear flows, 
the CRW approach cannot be used in the same way be- 
cause the basic building-block structures of CRW the- 
ory, namely kernel Rossby waves (KRWs)2, do not de- 
scribe buoyancy. However a formulation in the spirit 
of CRWs has been recently developed to include the ef- 
fects of stratification^. The basic building blocks of this 
generalized wave kernel approach, namely kernel Rossby- 
gravity waves (KRGWs), turn into KRWs when buoy- 
ancy effects are absent. 

We apply this approach to an example stratified shear 
flow in order to observe how the mechanics of this in- 
terpretive tool unfolds. The geometry is an infinite at- 
mosphere composed of divergence free fluids in a glob- 
ally constant shear flow. The atmosphere is comprised 
of fluids of three densities where the two density jumps 
of equal magnitude are located symmetrically away from 
the (nominal) midplane of the atmosphere. 



The important matter to recognize is that because 
there is no-basic vorticity gradient there will be no 
Rossby waves present at the density interfaces. In this 
way the KRGW building blocks of the generalized theory 
will reduce to what we call here kernel gravity waves or 
"KGWs" for short. The KGW concept has antecedents 
tracing back to the work of Sakai (1989) 7 in which 
density-shear disturbances in a shallow-water model are 
analyzed in terms of physical wave coordinates. The 
problem investigated by Sakaii cast the emerging insta- 
bilities in terms of a wave interaction mechanism in which 
counter-propagating waves influence each other from a 
distance. Baines and Mitsudera (1994)£ further gener- 
alize these concepts as fitting into a more unified mech- 
anism of shear flow instabilities. As an example they 
show how the classic Holmboe instabilityi^ may be ra- 
tionalized in terms of the interaction between the waves 
on separated density and vorticity surfaces. 

In the problem we examine here there will be four 
KGWs, two associated at each interface. Of the two at 
each interface one will travel faster than the local flow 
speed while the other will be slower than the local flow 
speed. Of the four modes only two modes, each counter- 
propagating with respect to the local flow speed of its 
respective interface, may transit into a stable-unstable 
pair. We find that for given values of the bulk- Richardson 
number (defined below) there will always be a band of 
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horizontal wavenumbers which admit normal-mode insta- 
bility. The KGW analysis performed here shows some- 
thing quite interesting: the mechanics of the instability 
bears strong resemblance to the mechanics responsible 
for the instability in the classic CRW problem^ and, fur- 
thermore, become mathematically equivalent in the limit 
of small horizontal wavelengths. Thus despite the fact 
that in this problem there are no Rossby waves present 
(i.e. as edge waves) buoyancy oscillations in relative mo- 
tion with each other interact dynamically in such a way 
that it can induce instability in the same way that CRWs 
produce instability in Rayleigh's classic problem 3 -^. 

Further reflection also shows that onset of instabil- 
ity satisfies the criterion of Hayashi & Young (1987)^. 
The stability criterion states that for flows of this sort 
a transition into linear stability occurs if there are indi- 
vidual waves in the flow which: (i) propagate opposite 
to each other, (ii) have almost the same Doppler-shifted 
frequency and (iii) can interact with one another—. In- 
deed in the current problem we see that instability de- 
velops out of KGWs which counter-propagate at their re- 
spective interface and through their action-at-a-distance 
interaction— trigger a transition into instability for given 
values of a bulk Richardson number. 

We note that the work of Caulfield 9 -, which is an in- 
vestigation of a variation of Holmboe's problem^ -, con- 
tains as a special limiting case the instability uncovered 
here. In that study a density configuration exactly like 
the one considered here is analyzed. The difference is 
that there exists symmetrically placed jumps in the back- 
ground vorticity resembling the classic Rayleigh profile. 
These vorticity jumps, located at z = ±d, are at points 
some distance removed from the density jumps located 
at z = ±ft. In general many types of modes appear in 
this problem. Of these, the ones referred to as Taylor 
modes may be shown to limit to the type of modes stud- 
ied here when d 3> h. However, this limiting procedure 
was not performed and the qualitative significance of this 
instability had not been appreciated at the time. 

It is important to note here that the three-layer den- 
sity configuration we consider in this work has been pro- 
posed to be qualitatively applicable to the question of 
drag reduction in rough seas driven by hurricane force 
winds 13,14 . The middle "foamy" layer— serves to dis- 
sipate energy and transfer momentum between the air 
and sea. Theoretical work in problems with this sort 
of stratificatio n 15 ! 16 consider uniform velocity profiles 
within each of the three layers but otherwise different 
between them. In this way we see the physical results 
derived from this study as complementary to these cited 
works because the velocity profile we consider has uni- 
form shear as opposed to (effective) delta-function jumps 
in the shear. 

This work is organized according to the following 
programme. In Section II we present and review the for- 
mulation derived in Harnik et al. (2007) 6 . In Section III 



we work through the theory for KGWs in an atmosphere 
of constant shear and a single density jump. As a matter 
of review we expend some effort describing the mechanics 
of the KGW. In Section IV we introduce the problem of 
two density jumps of equal magnitude, analyze its nor- 
mal modes, motivate its generalized formulation as an 
initial value problem and then argue for the rationalizion 
of its unstable dynamics. We find that when the system 
is unstable the mechanics of the modes closely parallels 
the processes occurring for unstable CRWs 3 under con- 
ditions where the layers are not too close to one another. 
We also spend a few words to rationalize the instability 
in terms of the criterion of Hayashi & Young (1987) and 
Sakai (1989)^11. In Section V we study the normal-mode 
response of this configuration when the layer is consid- 
ered foamy. This corresponds to the problem considered 
in Section IV but where, instead, the density jumps are 
not of equal mag nitude. As in Shtemler et al. (2007) 15 , 
the density ratio of successive layers is measured by the 
parameter S in which < S < 1. In comparison to the 
results of Section IV (where 5 w 1), we find that the 
range of unstable wavenumbers and the peak growth rate 
shrinks as 6 approaches zero. We also find that the modes 
are propagatory when they become unstable. Section VI 
summarizes our results and we show that the instabil- 
ity does not violate the Miles-Howard Theore m 19 ' 20 . We 
also suggest that this instability, manifesting itself un- 
der conditions that are classically considered to be sta- 
bly stratified, may be the simplest one possible which can 
destabilize a plane- Couette profile. We conjecture upon 
the results of this study and its relationship to the prob- 
lem of sea surface foam layers observed in hurricanes. 



II. BASIC FORMULATION: A REVIEW 

We begin the analysis by recasting the equations of 
linearized motion in terms of the formalism outlined in 
Harnik et al. (2007) 6 . Namely, the primitive equations 
of motion describing Boussinessq incompressible 2D flow 



Du 
~Dt 
Dw 
~Dt 



-U z w 

1 dP 

Po dz 



1 dP 

Po dx ' 



(1) 
(2) 



together with the equation of continuity and incompress- 
ibility 



Db 
~Dt 

du dw 
dx dz 



-wN 2 , 



0. 



(3) 
(4) 



The variables u and w are the horizontal (x) and vertical 
velocities (z) (respectively). The horizontal mean flow is 
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U and where we have defined U z 



as the mean shear 



together with -Sr = (d t + Ud x ). The density p is cast in 
terms of the variable b = —gp/po and where 



N 2 = - 



g_dp 

Po dz 



db 

dz 



Because the flow is incompressible the evolution of the 
density p(z,t) may be equivalcntly traced by the motion 
of some material invariant. In this case we follow the 
vertical displacement £ of a fluid parcel which was other- 
wise at rest. The undisturbed density profile is assumed 
to have a vertical variation of some sort and, thus, the 
displacement may be easily associated with clearly iden- 
tifiable surfaces. Thus, as was demonstrated in Harnik 
et al. (2007)£, we shall write the equations of motion in 
terms of the evolution of the vorticity and vertical dis- 
placement as 



Dq 
Dt 
£S 
Dt 



db 

-q z w+ — , 
ox 



( \ 

— = w{x,z), 



where q z 



U zz , b = —b z C, and where, 

dip , x dip 



w(x,z) = —, u(x,z) 



dz ' 



(5) 
(6) 

(7) 



in which we have expressed velocities in terms of the usual 
stream function formulation of 2D incompressible prob- 
lems and where the vorticity q is 



dw du 
dx dz 



v 2 ip 



(8) 



III. SINGLE DENSITY JUMP 



From here on out we assume the background flow state 
to be a constant shear (i.e. q z = U zz — 0). We examine 
here a single density interface located at the position z — 
0. The basic state density is written as 



p = po + ApH(z) 



(9) 



where H is the Heaviside function. The equations of 
motion now become, 

(d t + Ud x )q = {gAp/p )d x (6(z), (10) 

(d t + ud x )( = w(x,o). (ii) 

(flOl) says that the vorticity ought to also have a delta 
function dependence. Thus we make the ansatz 

q = q(x,t)5(z). 

To proceed, we assume that all variables have a Fourier 
decomposition in the x direction, that is to say 

~ e ikx + c.c. 



where k is the horizontal wavenumber. We begin with the 
solution to the streamfunction. As the ansatz indicates, 
away from the interfaces the vorticity is zero. Thus it 
means that one formally writes the solution to the stream 
function equation ([5]) as 



%j) = q(x,t) / G(z,z')dz', 



(12) 



since V 2 G = (d z — k 2 )G{z, z') = S(z — z ), and because 
this geometry is so simple the solution for tp is 



C L e -fr|*-*oi 

2k 



(13) 



Implicit in the construction of this solution is that (a) 
the vertical velocity, w = d x ip is continuous across the 
interface and (b) the vorticity delta function results from 
the jump in u across zq. We are now in a position to put 
these expression explictly into the equations of motion 
above. Thus one has, after shaking out algebraic factors 
and the solution ansatz, 



{dt 
(dt 



ikU )q 

iku )C 



Po 



(14) 
(15) 



where Uq = U(z = 0). Let us suppose that Uq = 
for this problem. Let us define the gravity wave speed 

c g = N/k = y^Hf • Provided that Ap < (stably 
stratified) then the gravity wave speed is real - and this 
we assume henceforth. We find that the above two equa- 
tions appear rewritten as 



d t q 
dtC 



-2ik 2 c% 



(16) 
(17) 



The above is the ode system, in time, describing the evo- 
lution of simple classical Rayleigh- Taylor modes. Note 
that this system is normal as far as the operators are 
concerned^. We can further proceed by defining new 
quantities through the relationship 



revealing to us the simplified system 



{d t + ikc g )C+ = 0, 
{d t -ikc g )C- = 0. 



(18) 



(19) 
(20) 



A normal mode analsysis may proceed from here by as- 
suming e~' lkct dependence upon all the solutions in which 
c is the wavespeed. It is quite straightforward to estab- 
lish that the eigenmodes for this system are described 
by 



{C-Cg)(+ 

(c + c ff )C- 



0, 
0. 



(21) 
(22) 
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For the sake of this illustration let us consider the propa- 
gation of the C+ mode which propagates at speed c = c g . 
To study C+ in isolation means setting f_ to zero which, 
in turn, means that, 

In other words, the disturbance height is positively cor- 
related with the local vorticity. 

Referring to Fig. [T]we depict this situation for a stably 
stratified density profile. We focus our attention to the 
section of fluid located between the two triangles (rep- 
resenting the peak and trough of the wave). The fluid 
layer has been disturbed from its equilibrium configura- 
tion (denoted by dashed lines) in such a way that it looks 
like a lever arm out of balance with the fulcrum denoted 
with a solid filled circle. One's intuition instructs that in 
this condition the lever arm has the tendency to rotate 
in a counter clockwise sense (positive vorticity) about 
this fulcrum point. This is simply the fluid's tendency 
to flatten out. Thus the out-of-equilibrium arrangement 
equates to a local creation of vorticity, i.e. a positive time 
rate of change of vorticity there (shown with a dashed 
counter-clockwise arrow on the graph). However, the 
already preexisting vorticity of this disturbance (shown 
with solid counter-clockwise arrow) causes, for example, 
the position of the fulcrum point to move up because the 
velocity field there points upwards (solid arrow). As this 
section of fluid responds, vertical velocity fields also de- 
velop at the positions of the triangles in proportion to the 
amount of vorticity being generated at the fulcrum (de- 
noted with the dashed vertical arrows above and below 
the left and right triangles respectively). 

In the panel immediately below we see what happens 
after one quarter cycle: the positive time rate of change 
of vorticity at the point which was once the fulcrum has 
turned into a positive vorticity and the fulcrum point has 
turned into the maximum height of the wave. As this 
occurred, the position of the triangles have moved corre- 
spondingly. In this way the entire pattern has shifted to 
the right. The mechanics of the oppositely propagating 
mode, C- 1 are recovered in a simple way: q and cor- 
respondingly w are reversed and hence dtq and dtC are 
then also reversed and then, consequently, the pattern 
moves to the left. Note that in this case (_ and q are in 
phase. 



IV. TWO DENSITY JUMPS IN A GLOBALLY 
CONSTANT SHEAR 

We now consider a situation in which there are two 
density jumps located symmetrically at a distance h away 
from the position z = 0. Thus, 




t = t D + Yd period 



pattern shift 
rightward 




FIG. 1: A mechanism schematic for the single kernel gravity 
wave in which positive disturbances are correlated with 
positive vorticity disturbances. Fluid of two differing densi- 
ties are disturbed with sinusoidal profile. The correspond- 
ing negative/positive vorticity profile is represented by solid 
clockwise/counterclockwise arrows. Restricting attention to 
between the two labeled triangles in the top panel, the density 
arrangement looks like a displaced lever arm in a gravitational 
field whose ID pivot (fulcrum) point is located at the position 
of the solid circle. This physically represents the tendency for 
the fluid to flatten out. The lever arm effect has a tendency 
to create positive vorticity about the fulcrum (dashed coun- 
terclockwise arrow). Concurrently, the preexisting flow due 
to the disturbance profile (solid clockwise/counterclockwise 
arrows) causes this point to rise upwards (solid arrow). As 
vorticity is created there, it also begins to create a velocity 
field at the triangles causing those positions to move (accord- 
ing to the manner denoted by the dashed arrows). By the 
time the fulcrum point reaches its peak, due to the lever arm 
effect there now exists positive vorticity in this local patch of 
the fluid after one quarter period. The pattern consequently 
shifts rightward (bottom panel). 



where the locations zi j2 are given by 

z 1 = h + d(x,t), z 2 = -h + ( 2 (x,t). 

We assume a configuration in which the shear is globally 
constant so that U = Az implying that q z — 0. We 
further say that A = U/h where U is the magnitude 
of the background flow at z = ±h. We have written 
matters in this way to always make sure that at z = h 
the wind is easterly (positive), in other words, we will 
assume that U > without loss of generality. Also, as 
before, we shall consider stably stratified configurations 
so that both Api < and Ap 2 < 1. 

Inspection shows that this form of the unperturbed 
density spawns delta function forms for the vorticity as 
we encountered before. In other words it means that we 
write 



p = ApiH(z - zi) + p + Ap 2 H(z - z 2 ) 



(23) 



qi(x, t)8(z - z{) + q 2 (x, t)8(z - z 2 ) 



(24) 
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Putting these forms into the governing equations ([5][6]) 
and assuming the single-wave Fourier ansatz 

{qt , Ci , 92 , C2 , ^} T = { 9i , Ci , 92 , C2 , 4>} T e lkx + c.c. , (25) 
we find the following four equations to emerge, 

Ft +ik 



ikhNf Ci, 


(26) 


w(h), 


(27) 


ikhN^2, 


(28) 


w(-h). 


(29) 



9 h P ^< 2 . We refer the reader to Appendix A 



with N 2 

1, 

for a detailed derivation of these equations from first prin- 
ciples. The solution to the streamfunction, i.e. (|A10|) . is 
rewritten here: 



1l k\z-h\ _ 12 k\z+h\ 

2k 2k 



(30) 



Since w — > ikip we may explicitly write the vertical ve- 
locity at the levels z = ±ft 



w 2 2 



Thus, we have finally 
d 



w{-h) 



2 2 ' 



jfi + ikUjii! = -2ik 2 c 2 gl Ci, (31) 
| + Hcr = -4e--4, (32) 



(9/; 



ikU C: 



-2ifc 2 c* 2 C2, 



..92 _^ e -2„ 

2 2 



(33) 
(34) 



in which 



_ gApi 



2p k 



"s2 



2po^ 



We may recast the equations of motion in this case 
in terms of the individual "kernel" waves that propa- 
gate on each interface (see Section III and Harnik et al. 
2007—). This re-formulation will be beneficial because 
it will show us transparently how these individual waves 
interact with one another. To that end we define the cor- 
related/anticorrelated ("±") wave quantities appropriate 
to each surface, 



C 2 ± = \ [ C2 ± 777—92 



2kc 



</2 



(35) 
(36) 



substituting these definitions into (|38|) along with the 
observation that 



9i 



2kc , 

si 



C 1+ -C-, 



92 

2fcc 



C2+ C2- ) 



and after some reshuffling of terms this results in, 



d 

— + ifecgi + ifcC/ ) Ci 
d 

— - ikcgi + ikU ) d_ 

— + i/cc ff2 - ikU ) C 2+ 



at 



- J, (37a) 

-»*A(C a+ -Ci_),(37b) 
-^A(Ci+-C 1 _),(37c) 
ifcc 92 - ifcf/J C 2 _ = -ife/3i(c i+ -C_),(37d) 



where 



-2fc/i 



n -2kh 



02 = C 



Provided that c g \ and c g2 are real (that is to say, that the 
density succesively gets smaller with increased height) , 
the structures £ 1± , £ 2± are each interpreted as pairs of 
kernel gravity wave modes at their respective surfaces. 
Note that £ 1+ and £ 2 _ are waves which move in the same 
direction as the local flow speed (at their respective lev- 
els) while Ci_ and £ 2+ move against the local flow speed. 
The latter two structures are central actors when insta- 
bility develops (see below). 



A. Normal Modes - and Recasting 

Now we can take this further and assume e~ lkct so- 
lutions and inquire about the normal mode behavior of 
(|37|) . revealing 

(c - c gl - U)( 1+ = (3 2 (4 + - C 2 _) , (38a) 

(c + c gl - COCx- = - C,-). (38b) 

(c-c g2 + U)( 2+ = j8i(c i+ -Ci_), (38c) 

(c + c g2 + U)( 2 _ = fa (C 1+ - Ci_) • (38d) 

One can quite easily generate a dispersion relation for 
the above set 

c 4 - (<& + c 2 g2 + 2U 2 )c 2 + 2U(c 2 g2 - c 2 gl )c + 

£/ 4 - U 2 (c 2 gl + c 2 g2 ) + c 2 gl c 2 g2 (1 - e~ ikh ) = 0. (39) 

From here on out we shall consider the simple profile 
c g i = c 2 2 = c 2 ., i.e. that the drop in the density profile 
is uniform from one step to the next, i.e. that Api = 
Ap2 = — Ap where Ap > 0. We write the gravity wave 
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speed as c 2 which is now assumed to be greater than zero 
(i.e. a real wavespeed). The solution to the wavespeeds 
becomes 



(l±C) 2 ±C 2 



id' 




(40) 



where K = kh and we have defined the scaled gravity 
wave speed with respect to the local shear velocity as 



C 2 = 



_9- 

u 2 



Ri 



Ri 



( 9hAp \ 



(41) 



where Ri is interpreted as the (bulk) Richardson number. 
We defined an "induced" wavespeed Ci n( j and its signif- 
icance will be explained shortly. In general this system 
has four modes given by ±c ± which are depicted in Fig- 
ure [2] for bulk Richardson number Ri — I. Inspection 
of the above shows that two roots are always real (the 
"+" branch) and we identify these as being propropagat- 
ing modes while the negative branch can admit an ex- 
ponentially decaying/growing pair and we refer to these 
as counterpropagating modes. This exists for a band of 
gravitywave speeds i.e. that is when the condition 



is satisfied 



1+C Z - 2Cyi 

and is so when 
K 



<0, 



K 



l-e- 2 « >Rl> TT^- (42) 
Figure [3] diagramatically encapsulates the result in ex- 
pression (|4"2"|) . There are two things to note. First is that 
for any given Ri there always exists a band of wavenum- 
bers supporting an instability. Secondly, inspection of 
(j42l) quickly shows that this band centers around the 
line Ri = K and that, further, the size of this band be- 
comes exponentially small as K ^> 1. In other words, 
as K — * oo the range of instability occurs within a 
small region centered on K whose width is measured by 
6 = e- 2K , i.e. 



K(l -6) <Ri< K{1 + 5). 

We ask the question: for the unstable mode in isola- 
tion, what are the relative amplitudes and phases of the 
four kernel waves with respect to each other? Let us say 
define the vector V = , ( 1 _ , £ 2+ , ( 2 _ ) T , then we see 
that (1551) can be rewritten as 



cV = LV 



in which /3 = c g e 2 /2 and 



Cg + U 
I -Cg + U 

P ~P 

P -p 



P 
P 



c g -U 



-p 
-P 



-Cg-U 



(43) 



(44) 



<B 1 



propropagating 



z counterpropagating 



" counterpropagating 



propropagating 



counterpropagating 



1 1.2 
K 



0.2 
0.1 

I 
-0.1 

-0.2 




FIG. 2: The dispersion relationship (|39[) f° r the case where 
c g i = c g 2 with Ri = 1, i.e. (|40|l . The top graph shows 
the real wavespeeds of the four modes as a function of K. 
The propropagating modes always remain stable while the 
two counter-propagating modes become unstable when they 
phase lock with each other and their phase speeds become 
zero. The bottom graph shows the imaginary wavespeeds of 
the counter-propagating modes. 



In the diagonal basis (denoted by primes) this system 
may be rewritten as 



cV' = CV, 



= diag(- 



,+c.,-c + ). (45) 



where V = WV for the rotation vector R. The columns 
of the inverse rotation matrix R -1 are composed of the 
eigenvectors of L (respecting the order of the eigenvalues 
appearing in C). We have written the matrix C so that 
the first position represents the mode that can go unsta- 
ble. Thus, if we wish to observe this mode in isolation 
we must require the following three equations 



0. 



0, 



hnVn = 0, (46) 



(n = 1,2,3,4) be simultaneously satisfied. (Note that 
this is a generalization of the wave-mode isolation pro- 
cedure we performed in Section III to follow the single 
kernel wave £+■) Since these are three equations with 
four unknowns, we ask what are the relative amplitudes 
and phases of the remaining three kernel waves when one 
has been set fixed to 1 . In this particular case of the un- 
stable mode, we find that the relative behavior is most 
transparently appreciated when we gauge the action of 
the kernel waves with respect to the kernel wave - 
i.e. the mode which, in the limit of large separation, cor- 
responds to the wave propagating against the mean flow 
at the upper interface, cf. ([38b ) with kh ^> 1. By choos- 
ing £j_ = 1 we can solve (|46|) for the remaining waves. 
Although M _1 may be written down analytically the so- 
lution for R is quite cumbersome and, as such, we depict 
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kh 



FIG. 3: The incidence of linear instability (shaded region) 
for given values of the bulk-Richardson Number Ri and the 
wavenumber k. In this particular example c ± = c g2 . Shown 
with a dashed horizontal line is the specific value of Ri = 1 
that is examined in the text. Instability occurs approximately 
for the range 0.80 < kh < 1.11 (shown with dotted lines). The 
range of instability shrinks as kh 3> 1 and stradles the line Ri 
= kh (also depicted). 

a typical example numerically of the resulting behavior 
in Fig. [3J To be specific, we write each kernel wave in 
terms of its amplitude and phase, i.e. 

By construction A 1 _ = 1, e 1 _ = 0. In the figure we 
plot the relative amplitudes and phases for the remain- 
ing waves with respect to £i_ as one sweeps through 
wavenumber K. 

We first note that in the limit of large wavenumber 
f3 becomes negligible and the four equations of ([T5|) de- 
couple. This means to say that the four kernel modes 
propagate not "knowing" of the other waves or inter- 
faces. This trend is represented in Fig|4] for kh ^> 1 
- which shows the amplitudes of the remaining modes 
(i.e. Ci+: C 2 ±) to go to zero as well. In the unstable zone 
(roughly 0.8 < K < 1.2) the amplitudes of ( ± _ and ( 2+ 
stay equal, i.e. A 1 _ — A. 2+ = 1, and the same goes for 
the the other two modes A 1+ — A 2 _ although they do 
not remain constant in the instability range as do the 
other pair. For this particular set of parameters the am- 
plitudes of A 1+ , A 2 _ are quite small compared to A 1 _ and 
this feature can aid us in developing a mechanical ratio- 
nalization of the instability (see next section). We note 
that the phases also show an eastward tilt with respect 
to each other in the unstable regime: e 2+ begins to show 
a positive phase with respect to e 1 _ = as well as e 
showing a positive phase with respect to e 2 _ = — it (the 
latter, incidentally, stays fixed throughout). 




0.5 1 1.5 

wavenumber K 



FIG. 4: Phase and amplitude diagram for the unstable mode 
(+c_) in which Ri — 1. All kernel wave phases and ampli- 
tudes are measured against the wave fj_ which is set with 
A 1 _ — 1, e 1 _ = 0. The panels on the right shows the complex 
wave speeds for the unstable mode (+c_ ) as a function of K. 
The band of instability in K is delineated with dotted lines. 
For K S> 1 all jnode amplitudes go to zero leaving behind 
only the wave In the unstable regime £ 2+ shares the 

same amplitude as C, x _ but begins to show an eastward tilt 
(in the same sense of the shear) of its phase, e 2+ > as one 
sweeps thru K. In the unstable phase, the amplitudes of the 
other two modes (i.e. Ci+>C2_) ^ so lock together and there 
is also an eastward tilt of their phases with respect to each 
other in this range (i.e. of e 1+ with respect to e 2 _). Note 
that the amplitudes of A 2 - and A 1+ in this regime are small 
compared to the other two amplitudes. The circles represent 
the function /3sin(e 2+ — £i_). 

Due to the symmetries of the normal modes we note 
that the action and relative configuration of the kernel 
waves for the stable mode (i.e. —c_ ) is also represented 
in the figure with two modifications: the first of these is to 
set c_ — * — c_ while the second would be to interchange 
the roles between C 1 _ and £ 2+ and between £ 2 _ and 
In the exponentially decaying regime £ 2+ starts to show 
a westward drift of its phase with respect to . 

B. Formulation as an initial value problem 



The equations describing the general linearized re- 
sponse also offers insights which will aid us in rationaliz- 
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ing the instability uncovered here. Thus we consider the 
general temporal response of the equations of motion (|37p 
by adopting the strategy of Heifetz & Methven (2005) 3 
and inserting into those equations the assumed forms 



where we have exploited the fact that Z 2 



The 



C 2+ = Z 2+ (t)e^ + W , C 2 _=Z 2 _ (i)e fc »- W . (47) 



After equating real and imaginary parts we have the eight 
ode's: four for the amplitudes, 

Z 1+ = k(3 [Z 2+ sin(e 2+ - e 1+ ) - Z 2 _ sin(e 2 _ - e 1+ }}8&) 
K- = k P [ Z 2+ sin(e 2+ - £ ± _) - Z 2 _ sin(e 2 _ - e^jj^b) 



Z 2+ = k/3 [-Z 1+ sin(e 2+ - e u 



sin e. 



-Z 1+ sin(e 2 



Z 2 _ = k/3 
and four for the respective phases, 
e 1+ = -k(U + c g ) 

r; 

-k/3 



Z 1 _ sin(e 2 



cos e, 



cos e, 



= -k{U~Cg) 



-k(3 



— cos(e 2+ - e x _ 



i 2+ = -k(-U + c g ) 



-k(3 



cos 



e 2 _ = —k(—U — c g ) 
Z^ , 



-k/3 



■ COS £, 



- 



— cos(e 2 

— — COS £, 



— cos(e 2 



- e, 



(48e) 



(48f) 



(48g 



. (48h) 



From the analysis of the previous section, an unstable 
normal mode implies a number of relationships to stand 
between the quantities appearing in (|48[) : that (i) that 



= 0,e 2 



since the modes do not "propagate" 



but are standing in that instance, (ii) e 2+ = 7r + e 1+ , and 
e 1 _ = 7r + e 2 _ , (iii) Z 2+ = Z 1 _ and Z 1+ = Z 2 _ . If wc 
assign to a the meaning of a growth rate then it is related 
to the phase differences according to 



a = k(3 sinl e 2+ — e 1 _ 



in Figure 1 we have depicted a - in terms of Im(c) which 
is equal to —a/ik. 

In the extreme limit where f3 <C 1, together with the 
above stated relationships, we may establish the relative 
scaling behavior between Z t _ and Z 1+ and between Z 2+ 
and Z 2 _ in the unstable normal mode regime. Taking for 
instance (I48c[) , if one is in the instability range then it 
immediately follows that 







U - 



o (p 2 ) , 



(49) 



scaling on Z 2 _ immediately follows from the equality of 
it to Z 1+ . We learn from this analysis an important clue: 
namely that in the limit where /3 <?C 1 the main actors 
involved in the development of the normal-mode insta- 
bility are the kernel waves C 2+ an d Ci- while the other 
two kernel waves, £ 2 _ and £ 1+ , should be viewed as being 
slaved (in a qualitative sense) to the former two. This 
fact allows us to interpret the instability with greater 
transparency. We devote our attention to this in the 
next section. 



C. The normalmode instability rationalized 

Although for this problem there are four KGW waves, 
the general normal-mode analysis indicates that expo- 
nential behavior is mainly rooted in the two individual 
KGWs which move against the prevailing (background) 
flow. Take for example the KGW propagating against 
the flow associated with the top interface: if one were to 
move into the rest frame at that level, i.e. in a frame 
moving eastward with velocity U, then this wave would 
appear to move westward. By symmetry reasons the cor- 
responding KGW on the bottom interface would appear 
to move eastward. Thus to an observer in the labora- 
tory frame the culprit KGWs propagate with a speed 
which is slower than the local flow. As mentioned ear- 
lier, we shall refer to these waves with this character as 
counterpropagating-KGWs. Taken as individual waves, 
i.e. free of cross-layer interactions, the counterpopa- 
gating KGWs are those KGWs whose wavespeeds are 
±(U — c g ). Those KGWs that propagate with the local 
flow, i.e. those which have wavespeeds equal to ±(C/+c 9 ), 
will be referred to as propropagating-KGWs. 

We further remind ourselves that when layers interact, 
the primary (eigen)modes leading to instability come in 
pairs - showing either an eastward or westward propaga- 
tion of the eigenmode pattern. In general the individual 
eigenmodes are structures composed of all four KGWs 
however we see from the arguments in the last section 
that when (3 is small, the eigenmodes are primarily con- 
structed of only the two counter-propagating-KGWs d_ 
and ( 2+ while the propropagating-KGWs, namely ( 1+ 
and £ 2 _ are (in a sense) slaved to the former two. In- 
stability develops by which the modes in question show 
wavespeeds passing through zero. Thus at marginality 
we have two wave patterns that both appear standing to 
the observer in the laboratory frame. 

From these observations we shall argue here that the 
mechanism for instability in this problem is strongly anal- 
ogous to the process leading to the instability of CRWs as 
in Rayleigh's classic problem. Thus we propose that the 
dynamics are governed only by the action of the counter- 
propagating-KGWs for <C 1 (i.e. K ^> 1). Since in 
that case the propropagating-KGWs are slaved to the 
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counter-propagating ones, we posit that the dynamical 
evolution of this system (for c g \ = c g 2 = c g ) is given by 
the more simpler form 



+ ikcj d_ = -ik/3( 2+ , 
- ike ) 4+ = ifc/3Ci- 1 



(50) 
(51) 



where c = U — c g . We note immediately that mathemat- 
ically speaking, (|50ti51|) have the same structure as the 
equations describing the evolution of interacting kernel 
Rossby waves (KRWs, cf. 9a-c of Heifetz & Methven, 
2005 3 ). It means to say that these mutual interaction of 
KGW waves (in this limit) may be considered and ratio- 
nalized in the same manner in which interacting KRWs 
are rationalized. The familiar instability properties ex- 
hibited in the problem of CRWs follows onto the quality 
of these dynamics as well, in particular, the co-action of 
hindcring-hclping of individual waves when instability is 
present. 

We may further rationalize this instability, as viewed 
through through the lens of the reduced model (150M51|) . 
in terms of the stability criterion of Hayashi & Young 
(1987) and Sakai (1989)L±i referred to in the Introduc- 
tion. Inspection of the dispersion relationship for the sys- 
tem (|50H51|) . e.g. the counter-propagating modes in Fig. 
EJ shows that the those modes which become unstable 
propagate opposite to each other for wavenumbers just 
beyond the onset of instability. Their Doppler-shiftcd 
frequencies at onset (i.e. when c = for both waves) 
are ±kc respectively and they are almost equal by virtue 
of c being almost zero in the limit appropriate for this 
reduced model. And finally without the term f3, which 
represents the mutual interaction of these waves, there is 
no chance for instability to occur. 



V. NORMAL MODES OF A FOAMY LAYER 

We consider further results assuming the two density 
jumps represent the configuration of a foamy layer as dis- 
cussed in the Introduction. In particular we will examine 
the atmosphere with stratification considered in Shtem- 
lcr et al. (2007). 15 In other words we shall analyze an 
atmosphere in which the density of the air p a is 5 times 
the density of the foam pf which is, in turn, 5 times the 
density of the sea layer p s . In terms of the density formu- 
lation given in (|23j) we make the following replacements 
po -> p f and 



Apr 



1 



Pf 



Ap 2 ^- Pf (l-S), 



and for the velocities we have 

,,2 



zL 
5 ' 



c 



</2 



2 _ 5(1-5) 

c f - 



where < 5 < 1 . Considerations of this system in a small 
vicinity near 5=1 recovers the dynamical behavior of 
the case studied in the previous section. The dispersion 
relationship (|39p now appears as 



- 2 + C 



2 l + 5 



1 + 5 C) 
f ~5 



c \ 2 o 1 — 5 ( c 

u) - 2C f—(u 



l-C?^ + ^(l-e- 4A ') = 0. 



(52) 



where Cf = Cf/U. Inspection of the above relationship 
suggests that when instability occurs, the unstable waves 
will have some propagation associated with them unlike 
the stratification considered in the previous section. In 
a similar fashion we define the quantity Rf such that 
Rf/K and where 



c 2 



l 



R f = 



gh{l-5) 
2U 2 ' 



2k 



which is related to the inverse of the Froude number. 
We show in Fig [5] the general trends for these circum- 
stances when Rf = 1. A non-zero value of 5 means that 
when instability sets in, it is no longer stationary as in 
the previous section, but the modes in question are now 
propagatory. In the range of unstable wavenumbers the 
propagation speed is is also a function of K. For Rf 
held fixed, as 5 approaches zero the range of unstable 
wavenumbers gets smaller, the location of the unstable 
band shifts to larger wavenumbers and the peak growth 
rate in this unstable band similarly get smaller. 

In Fig. [5] we display the results of fixing 5 and varying 
the bulk Richardson number Rf. This time, however, 
we use numbers that figure to be relevant to the prob- 
lem of hurricanes. Because the density ratio between 
water and air is one thousand to one we assume that 
5 = V0lX)T « 0.03. Following Shtemler et al. (2007)^ 
we take (as a rough order of magnitude estimate) for the 
foam layer thickness the figure h = lm. Because, empiri- 
cally speaking, strong winds over water follow a logarith- 
mic profile Powell et al. (2003)^ generate fits of this form 
to the data accumulated during hurricane over-flight mis- 
sions. For hurricane force winds (i.e. winds exceeding 34 
m-s -1 at 1km above the sea surface) the wind speed at 
lm above the sea surface is approximately 20 m-s -1 as 
inferred from Figure 2 of Powell et al. (2003)^4 Thus 
we assume this figure for U. Roughly speaking, it means 
that the bulk Richardson number for these conditions is 
Rf ~ 0.01. The results depicted in Fig. [6] are for three 
values of R f = 0.04,0.01,0.0024. The trends appearing 
are that the peak growth rate increases as Rf is made 
smaller although the unstable range in K begins to nar- 
row as well. Just as in the results of the previous section 
(e.g. see Fig. [3]) there exists a critical value of Rf(5) 
below which K = is also unstable. In general this is 
probably a pathology of the theory and one ought not to 
take too seriously the results at K = 0. However, once 
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FIG. 5: Wavespeeds and growth rates as a function of 
wavenumber K for various values of 8 at Ri = 1. The top row 
of panels demonstrates the wavespeeds. Whenever instability 
sets in the wavespeeds of two modes merge. Unlike S ~ 1 
theory, the wavespeeds are non-zero and have K dependence 
when instability sets in. The bottom panel shows how the 
growth rate behaves as S becomes small: (i) the range in K 
for which growth occurs narrows and, (ii) and the amplitude 
of the maximum growth rate reduces. 



R f dips below this value the range of unstable wavcnum- 
bers narrows. We note finally that the first wavenumber 
to be unstable at Rf — 0.01 occurs at K m 0.15. We 
shall make a remark about this and its relationship to 
hurricane data in the next section. 



VI. SUMMARY AND DISCUSSION 
A. Recapitulation 

We have studied how KGWs interact across two lay- 
ers in a medium with globally constant shear. KGWs, 
along with KRWs^, are dynamical phenomena forming a 
subclass of KRGWs developed in Harnik et al. (2007)2.. 
In terms of the formalism developed, we have found that 
the interaction of KGWs under these conditions bear sig- 
nificant similarity to KRWs especially when instability 
sets in. For KGWs, the presence of a jump in the back- 
ground density at some layer means that these places sup- 
port the creation of vorticity which ultimately relates to 
the interpretation that such configurations are sources of 
baroclinic torques 7 . Unlike KRWs, where a jump in the 
background shear produces a single propagating Rossby 
edge-wave, an isolated jump in density creates two sur- 



FIG. 6: Wavespeeds and growth rates for a foamy layer as 
a function of wavenumber K for various values of Rf at S = 
0.03. The top row of panels demonstrates the wavespeeds 
of those modes demonstrating growth and decay only. The 
peak growth rates get larger for smaller values of the bulk 
Richardson number. 



face gravity waves propagating with equal and opposite 
directions. These cohabitating (in the sense that they ex- 
ist on the same layer) counter-propagating gravity waves 
can be cast into a a set of new variables which explicitly 
shows that the wave pair do not interact. 

The KGW kernel is composed of the measure of the 
edge-wave vorticity and the amount which the surface 
has been displaced from equilibrium. Viewed in terms of 
this combination of physical variables, we have demon- 
strated how a single wave mechanically operates as it 
propagates. The vorticity field corresponds to a velocity 
field around the disturbed surface while the disturbed 
surface indicates how further vorticity is to be generated 
ahead of the peak of the wave. The disturbed surface 
creates this vorticity (in the sense of the mechanical de- 
scription we present) because it looks like a displaced 
lever arm - when let free such a lever arm has the ten- 
dency to rotate about its equilibrium point and generate 
vorticity as a consequence. The true wave executes these 
steps continuously in a fluid way. 

The analysis of the interaction of KGWs on two lay- 
ers was facilitated here by assuming the density jumps to 
be the same across each layer which results in a normal 
mode problem that is analytically tractable. Because of 
the shear profile, the top layer moves with speed U while 
the bottom layer moves with speed —U. In this case 
one has four normal modes (a co-habitating pair on each 
layer). Normal mode instability arises primarily through 



11 



the interaction of two modes from each layer - these 
modes are ones which, in the absence of mutual inter- 
action, would propagate against the background flow of 
the layer. The remaining two modes contribute less and 
less to the development of the instability when the layer 
separation (or, equivalently, the horizontal wavenumber) 
gets large. In this large separation limit, the resulting 
effective equations have the same mathematical form as 
the ones describing the (exact) development of unstable 
CRWs 2 -. This means to say that the mechanism behind 
this instability is qualitatively similar to the way insta- 
bility emerges in the problem of CRWs. 



B. As being a CRW analog 

In a sense the instability is the analog of the CRW 
instability and it is for this reason why we refer to this 
mode interaction process as the instability of counter- 
propagating gravity waves. Indeed the KGW may be 
interpreted as a (delta-function) vortex structure which 
propagates at some speed on the layer in question. 
Whereas in KRWs the source of the vorticity is the 
jump in the background vorticity, in KGWs the source is 
the baroclinicity that arises when the jump layer is per- 
turbed. The ensuing dynamics are qualitatively similar 
in all other major respects. It is also true that both left 
and right going waves become excited on a layer when a 
single KGW propagates along the other. In this way the 
conception and rationalization of the dynamics of four 
KGWs is more complicated. But, and we reiterate this 
point, the basic qualitative features of the instability in- 
volves the dynamical influence between only two of the 
KGWs, namely those that counterpropagate with respect 
to the background flow at their respective interfaces, and 
the processes at work resemble the mechanism responsi- 
ble for the instability of CRWs. 



C. The normal- mode instability, its rationalization 
and relationship to the Miles-Howard Theorem 

As the disturbances considered here preclude the pos- 
sibility of an inflexion point instability (i.e. the Rayleigh 
criterion), the character of disturbances falls under the 
purview of the Miles-Howard theorem. The Miles- 
Howard theore m 19 ' 20 for the stability of plane-parallel 
stratified shear flows states that a sufficient criterion 
for stability of normal-mode disturbances is that the 
Richardson number (J) of the flow be everywhere greater 
than j in the domain. That is to say that a necessary 
condition for instability is that somewhere in the flow 

T = _l^P 1 „ 1 /„, 

J ~ pdz ' {dU/dzf < 4' [M) 



where the quantities defining J are in dimensional units. 
Because this is a plane- Couette flow dll/dz is everywhere 
defined and not equal to zero. For the "stably" strati- 
fied fluid we have considered here, i.e. — (g/p)dp/dz>0, 
J would then be described by two delta- functions, with 
positive coefficients, centered on the locations of the den- 
sity discontinuity. This means to say that J is zero ev- 
erywhere except at the points where the density discon- 
tinuity occurs. At the very least, the instability uncov- 
ered here satisfies the necessary Miles-Howard criterion 
for normal-mode instability as J is less than j in sub- 
stantial parts of the flow. 

What is interesting to reflect upon is that this insta- 
bility might be the simplest kind possible in a shear flow 
absent vorticity waves. Regular (inviscid) plane-Couette 
flow, though satisfying the necessary condition for insta- 
bility, is however stable as is a plane-Couette flow with 
a single density interface. All that is needed for a plane- 
Couette profile to exhibit normal-mode instability is for 
there to exist a second density interface. 

Given the previous discussion it seems to us that 
normal-mode instability in shear flows merely requires 
the presence of counter-propagating waves of any sort. 
Indeed a cursory review of some previous work re- 
veals: (i) the classic instability emerging from the model 
problem investigated by Rayleigh^ has been reinter- 
preted as arising out of the interaction of phase-locked 
counter-propagating vorticity waves^ 2 - (i.e. Rossby edge 
waves), (ii) a simplified restricted version of Holmboe's 
instability 10 has been shown to be emerging out of the in- 
teraction of a pair of counter-propagating waves in which 
one is a vorticity wave and the other is a gravity wave® 
and, (iii) here we have demonstrated yet another insta- 
bility borne out of the interaction of counter-propagating 
waves which are comprised, in this instance, of phase- 
locked gravity waves. 

This leads to the apparently counterintuitive conclu- 
sion that a convectively stable stratification can destabi- 
lize a flow profile that is stable in the unstratificd case. 
It is nevertheless important to recall here that this result 
is not in contradiction with conventional linear stability 
criteria. Indeed, in the presence of stratification, one of 
these criteria is the Miles-Howard criterion found in |53|) 
which places no limitations on the wind curvature. It 
means that there is a much larger array of wind profiles 
that can be unstable in the stratified case than in the 
unstratificd one. 



D. Foamy layers and a conjecture 

We observe that the range of unstable wavenumbers is 
finite in the case where the density jumps are the same 
(Section IV). But, in comparison, both the wavenum- 
ber range and the maximum growth rate of the insta- 
bility when considered in a foamy layer (Section V) re- 
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duces as the the density parameter describing the foam 
8 gets small. However consideration of parameters hav- 
ing more relevance to real atmospheric phenomenon re- 
veals some interesting predictions. When numbers ap- 
propriate for observed hurricanes are utilyze d 14 ! 15 we find 
that they correspond to bulk Richardson numbers (Rf) 
in the vicinity of 0.01. Inspection of the growth rates 
show that the wavenumber for the onset of instability is 
kh w 0.15 in which h is the nominal height of the foam 
layer. Assuming, as we have, that h = lm it means that 
the wavelength of this unstable mode is approximately 
A ~ 45m. An inspection of Figure 4a in Powell et al. 
(2003) 14 shows a photograph of the sea from a quarter 
kilometer height during a hurricane whose windspeeds, as 
measured at the height of the aircraft, are approximately 
45 m-s . The image shows the foamy patches and the 
length scales of the groupings that appear is approxi- 
mately 30m in length. Figure 4b of the same work shows 
the foamy patches after the windspeeds have increased 
to 55 m-s -1 and one sees clearly that the foam patches 
have covered the entirety of the sea surface. Since larger 
values of the windspeed correspond to smaller values of 
critical wavenumbers, the theory developed here shows 
at least the same trend in the observed photographs. 

Of course there are obvious shortcomings of applying 
the theory developed here to the circumstances noted 
for this hurricane observation: (i), that the height of the 
foam layer was assumed since there are no direct observa- 
tions to that end from the data available, (ii) it is not at 
all clear if there exists a "laminar" foam layer from which 
a secondary instability develops creating the frothy foam 
layer observed in these photographs and, (iii) the assump- 
tion of a linear velocity profile in place of a logarithmic 
one is clearly a gross simplification. The latter assump- 
tion is probably satisfactory on the scale of the foam layer 
but certainly breaks down if one tries to extend it much 
beyond a logarithmic scale height of the layer thickness 
(e.g. 10m or higher if h = lm). Thus although the ap- 
parent consistency observed between this simple theory 
and the observations made is not proof that the theory 
developed here is relevant to such complicated processes 
(i.e. hurricane-sea interfaces), it is at least encouraging 
that this theory makes predictions which are of at least 
the same order of magnitude as the structures observed. 
As such we treat this discussion as being purely conjec- 
ture at this stage. 
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APPENDIX A: DERIVATION 

We show here the step-by-step derivation that ulti- 
mately leads to the equation set (|26H29[) . The basic den- 
sity profile and its derivatives, written in terms of 6, is 
itemized here. 
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these expressions are inserted into the governing equa- 
tions ([STo]) . 
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From inspection above it is clear that the vorticity of 
the system takes on delta function character at all of 
the jumps of the system. This then means applying the 
ansatz 

q = qi{x, t)6(z — z\) + q 2 (x, t)5(z - z 2 ). 

Collecting all terms of like delta functions, one may set 
each of the coefficients to zero because of the mutual 
"orthogonality" of the delta functions. In a sense this 
system has been discretized. This results in the following, 
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(A6) 
(A7) 
(A8) 
(A9) 



To recover (I26H29[) we apply the Fourier ansatz (f2"S"l) to 
(IA6IIA9I) . With which it immediately follows that the 
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streamfunction field ip (subject to the usual jump condi- 
tions) is the familiar form, 

i> = - ( %-e- k ^ + ^e- k ^ + 'A e lkx + c.c. (A10) 
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